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Abstract 

Projected momentum distributions of electrons, i.e. Compton profiles above the topmost atomic 
layer have recently become experimentally accessible by kinetic electron emission in grazing- 
incidence scattering of atoms at atomically flat single crystal metal surfaces. Sub-threshold emission 
by slow projectiles was shown to be sensitive to high-momentum components of the local Compton 
profile near the surface. We present a method to extract momentum distribution, Compton profiles, 
and Wigner and Husimi phase space distributions from ab-initio density-functional calculations of 
electronic structure. An application for such distributions to scattering experiments is discussed. 
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I. INTRODUCTION 



The local density of states above the topmost layer, the occupied states of which represent 
the electron density spilling out into vacuum, is an important quantity characterizing many 
properties of the surface and image states l[ 0], surface magnetism 0, 0], and chemical 
reactivity mediated by charge transfer Ig]. Yet, many of its properties are difficult to 
directly access experimentally because they are overshadowed by bulk contributions. A 
case in point is the projected momentum distribution, the Compton profile of the electronic 



structure. While bulk Compton profiles 



lave been mapped out in considerable detail by 



electron scattering [rj, fast ion scattering ^, and X-ray scattering j^, little is known about 
the surface Compton profile characterizing the anisotropic momentum distribution of the 
ground-state electronic structure reaching out into vacuum. 

Recent experiments for grazing-incidence scattering of low energy neutral rare-gas atoms 
at Al(lll) and Al(llO) single crystal surfaces provide first evidence of high-momentum 
components in the electronic selvage jl^. The projectile threshold velocity vth required for 
the kinetic emission of surface electrons by atoms in a head-on binary encounter collision 
{2vth + Ve)'^ /2 > (Ep + W) {Ep: Fermi energy, W: work function) was found to be below the 
standard value v^j^ for Ve = vp {vp: Fermi velocity). This sub-threshold kinetic emission 
indicates the presence of off-shell velocity (momentum) components well above the Fermi 
velocity Ve = vp. Obviously, the frequently employed free-electron gas approximation for 
conduction band electrons is not sufficiently accurate to account for sub-threshold kinetic 
electron emission (KE) |lo|. The experimental data can be taken as first evidence for the 
presence of momenta above the Fermi momentum qp in realistic momentum distributions 
above the surface taking many-particle effects (correlation) and the crystal potential into 
account. 

In the present communication we discuss the extraction of the momentum distribution 
and related quantities from ground-state density-functional calculations of the surface band 
structure. For estimating impact-parameter dependent emission probabilities in grazing- 
incidence surface scattering information beyond the projected momentum distribution, the 
Compton profile, is desirable: the quantum phase-space distribution or Wigner function 



12j . p {q,-r). It allows to determine the correlation between position above the surface 



and local momentum encoded in the electron density compatible with the Heisenberg un- 
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certainty principle. We also consider the Husimi distribution |l3^, p^{q,r), which represents 
a minimum- wavepacket average over . 

We present numerical results for the Al(lll) and Al(llO) single crystal surfaces using 
the ABINIT code Q| with soft pseudopotentials Applications to KE will be discussed. 
Atomic units are used throughout unless stated otherwise. 

II. MOMENTUM DISTRIBUTIONS AND DENSITY FUNCTIONAL THEORY 

State-of-the-art electronic structure calculations for the bulk and surface invoke density- 
'unctional theory (DFT for electronic structure programs employing DFT see, e.g. 
3 l^, 13 )• Ground-state energies as well as excitation spectra can be calculated 
within linear-response with sub-eV accuracy. Apart from the ubiquitous uncertainty with 
regard to the appropriate exchange-correlation potential, V^c, application of DFT to other 
classes of observables faces the conceptual difficulty that appropriate read-out functional 
are, a priori, not known unless the observable can be directly expressed in terms of the 
electronic density p(r). In the present contribution we inquire into a read-out functional for 
the Compton profile, i.e. the momentum distribution p{q) near surfaces. 

DFT employs the solution of a set of one-particle Kohn-Sham (KS) equations for the 
pseudo-wavefunctions ipi 

Heffirmr) = (-^A + V,ff[p{r)]^ = SiM^ (1) 

which have the property to add up to the exact electron density 

p{r-) = J2\Mrl\\ (2) 

provided the exact functional form of the exchange-correlation potential Vxc entering the 
effective potential K// is known. This, however, is not the case and a variety of approximate 
functionals suited for specific problems are in use. 

The quasicontinuum of bands n filled up to a maximum energy Ep is replaced by a 
discretized set of wavefunctions with band indices n and k. Then, for any band with band 
index n the wavefunction associated with momentum k in the first Brioullin zone solving 
the Schrodinger equation for the solid can be written by a sum over plane waves 

^n,drl = J2c^^dG)e^C^-'^>^ (3) 

G 



where G is any reciprocal lattice vector. Employing periodic boundary conditions, first- 
principle codes compute the electron density either for infinitely extended crystals (compu- 
tational box is equal to the crystal unit cell) or, at crystal surfaces, by calculating the ground 
state for an infinite number of thin crystal slabs ("supercell methods"). In either case, the 
coefficients ^ are returned as output of the calculation. The total electron density is given 
by, 

k 

where occ^^ ^ is the occupation of the (pseudo-) state with indices n and k and is the weight 
for the point k in reciprocal space resulting from symmetry considerations if calculations 
are restricted to the irreducible BriouUin zone. In our case, = const as we work on an 
equal-spaced grid extended over the whole BriouUin zone. 




A. Momentum distributions p{q) 

For a Schrodinger wavefunction (j){r) coordinate-space densities and momentum densities 
are related through the Fourier transform of the wavefunction, 

Piq) = \mf (5) 

with 

0(g) = (27r)-3/2 / d^re-''^^{r) . (6) 



Since the KS pseudo-wavefunctions ^{r) are, in general, not to be identified with true 
wavef unctions, the applicability of Eqs. El and El is, a priori, not obvious. Note that p{r) 
and p{q) are not Fourier transforms of each other. As an alternative to applying Eq. (HI 
to KS-wavefunctions, it might be possible to formulate the Hohenberg-Kohn theorem in 
momentum space with p{q) as fundamental quantity. In this case the KS equations would 
become integral equations and Vxc a non-local integral operator. To our knowledge, this 
avenue has not yet been explored. 

In the following, we postulate the applicability of Eqs. El and El for KS pseudo- 
wavefunctions. This can be justified for the homogeneous electron gas. Its extension to 
realistic DFT calculations for surfaces can be viewed as an analogue to the local-density 
approximation to the exchange- correlation potential. 
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FIG. 1: (Color online) Momentum distribution of an Al(lll) slab. The Fermi momentum for Al 
is about qp ~ 0.9 a.u. Higher momenta represent correlation and crystal structure effects. 

Accordingly, ip^ ^ is given in momentum representation by 

G 

= (27r)-^/2 J2 ^n,dG) S{q- ik + G)) . (7) 
G 

For the three-dimensional momentum distribution we find 

v-^ I ~ |2 

n,k 

n,k G,G' 

= (27r)-35^occ„,:«;,^5^|c„^,KG)|'5(g-(fc + G)). (8) 

n,k G 

Fig. ^ shows cuts through the momentum distribution of an Al(lll) single crystal slab 
at fixed values. As Al is well described by a nearly-free electron approximation, the 
momentum distribution closely resembles a sphere with radius qp ~ 0.9 a.u. but features 
also small but finite contributions at higher momenta due to correlation and the crystal 
structure potential. Furthermore, calculating p{q) for a crystal slab leads to an additional 
increase of higher momenta normal to the surface (g^) due to the breaking of translational 
symmetry at the surfaces. 
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B. Distance dependent momentum distributions 



If we perform the Fourier transform only for coordinates parallel to the surface we retain 
the ^-dependence (coordinate along the surface normal) of the position distribution but 
extract the momentum distribution along the coordinates in the surface plane p{q\\,z). It 
can be rewritten as a sum over plane-wave coefficients and delta functions, 

n,k 

= i2nr^ Yl ^^^n,m E ^nS(G)clk^G') e^(«^-^'> x (9) 

n,fc G,G' 

X (fell (fell 

Fig. 121 shows the distance dependent |-distributions for Al(lll) (left panels) and Al(llO) 
surfaces (right panel) averaged over the azimuthal angle (pq as a function of distance z from 
the topmost atomic layer on a logarithmic scale (color coding). The Fermi momentum is 
indicated by the black arrows. Differences between the distributions provide information on 
the face dependence of the momentum distributions and on the underlying surface potential. 
For example, the hump near the wavenumber q = 27r/as ~ 1.16 a.u. in the top left panel 
of Fig. El originates from the nearest neighbor distance in the Al(lll) surface ~ 5.4 a.u. 
In the case of the Al(llO) surface (right hand panels) the existence of "wide" and "narrow" 
surface channels becomes evident by the different extension of the distribution in and qy 
directions (bottom right). 



C. Wigner distributions p^{q,r) 

The Wigner distribution is the quantum mechanical analogue of the classical phase-space 
distribution p'^\q,r). It features remarkable resemblance to p'^\q,r) as it fulfills the same 
relations as p^' 

p(g) = j d^fp^{q,f) (10) 

p(f) = I d'qp'^{q,f). (11) 

p^ is real valued. An interpretation in terms of a probability distribution is, however, not 
possible as p^ is not positive definite. The latter is the immediate consequence of the 
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FIG. 2: (Color online) z-dependent momentum distributions as a function of the distance from 
the topmost atomic layer. Top left and right panels show distributions for Al(lll) and Al(llO) 
surfaces, respectively. Black dots indicate experimental estimates for local Fermi momenta (see 
chapter IlIII) . Bottom panels show — Qy distributions for both faces taken at z = 1 a.u. above 
the topmost atomic layer (dashed line in top panels). The white arrow in the bottom left panel 
indicates the wavenumber q = 2n/as ~ 1.16 related to the nearest-neighbor distance of Al(lll) 
surfaces. 

position-momentum uncertainty, can be calculated from the Fourier transform of the 
off-diagonal elements of the single-particle density matrix 

n,k 

and is given by 

p^(g, r) = ^ y (f - y\p\r + V) e'^'^U'y . (13) 

Assuming now again that the KS orbitals can be taken to represent the one-particle 
density matrix, the calculation of using DFT input is straightforward. It should be 
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FIG. 3: (Color online) Wigner function p {q,r) in front of an Al(lll) surface, qz = 0, z = 1 a.u. 
above a surface atom. Around the main part of the distribution, negative valued areas (purple 
contour, also in Fig. 0]) impede the interpretation as probability distribution. 

noted that the use of KS-orbitals assures that the exact position density results from the 
reduction of the Wigner function (Eq.llip. This observation suggests that the representation 
of in terms of KS orbitals and thus also of p{q) (Eq. [TUI) is a meaningful approximation. 
For a single-band wavefunction ^ we find 

G,G' 
G,G' 

« E ^nrkiGK/Gy^'^-^'>^6{2q- 2k-G-G'). (14) 

G,G' 

The total Wigner function is therefore 

p^{q,r) oc '^occ^ f:w^p^0,r) . (15) 

n,k 

The failure of the naive interpretation as a probability distribution can be seen from Fig. 
El p^ always features negative regions and strong oscillations in the regions of interest. 
Fig. m shows cuts through the Wigner function for various positions above the irreducible 
part of the surface unit cell (see inset) at ^2 = a.u. As Wigner functions are a natural 
starting point for quantum calculations of scattering processes, details of the distributions 
may account for features in electron emission experiments sensitive to details of the projectile 
trajectory along the surface. 
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FIG. 4: (Color online) Wigner function for an Al(lll) surface evaluated at 9 points of the irre- 
ducible part of the surface unit cell (see inset) and Qz = z = 1 a.u. 



D. Husimi distributions p^{q,' 



The Husimi distribution is a convoluted (or averaged) Wigner function and even closer 
linked to classical phase-space distribution than the Wigner function. Husimi suggested to 
average the Wigner function over a minimum uncertainty wavepacket. is defined by 



p^iq^r) = — I d'^x I d'^p {p, x) exp 



{p- g)^ 



a' 



(16) 



where a determines the width of the wavepacket in q and r directions, p^ is the closest 
analogue to classical phase-space distribution consistent with the Heisenberg uncertainty 
principle. In particular, p^ is positive definite and allows for a probability interpretation. 
Eq. El would be time consuming to evaluate directly as it would require the evaluation of 
p^ for all coordinates q and r and performing a 6-dimensional integration. This, however, 
can be circumvented by resorting to the equivalent expressions 
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[x — r) + iqx > d x 



and 



p^(g,f) oc J]occ„_^w^ Pn,kiP^^ 



(17) 



(18) 



n,k 



which reduces the numerical effort considerably, a in Eq. El again determines the width of 
the wavepacket in r*. The averaging process smooths out oscillations of p^ on the scale of 
the de Broglie wavelength. It renders p^{q,r) positive definite over the whole phase space 



and therefore allows for the desired interpretation as the probability of finding an electron 
with momentum q at position r. It has, however, the disadvantage of featuring unphysically 
high momentum components due to the infinitely long tail of the Gaussian distribution in 
Eqs. El and El Employing such a distribution in scattering calculations will therefore lead 
to a spontaneous escape of electrons with energies higher than the surface potential unless a 
suitable projection formalism is applied. Alternatively, reducing the width of the wavepacket 
in momentum space would also diminish the weight of the high-momentum components but 
would come at the price of an averaging over a large volume in f space (and therefore z) 
smoothing out the density gradient at the surface. The Husimi distribution can therefore 
be mainly used for an intuitive interpretation of quantum mechanical results in terms of the 
quasi-classical probability distribution. 



III. APPLICATION TO ABOVE SURFACE KINETIC ELECTRON EMISSION 



Outside the topmost atomic layer the electron density is quickly reduced. Electron wave- 
functions reaching farthest out of the surface will have large momentum in direction normal 
to the surface {Qz)- Therefore, smaller momenta parallel to the surface are observed than 
in the bulk. This behavior has recently been experimentally investigated by Winter et al. 



21[. They have scattered atomic projectiles off single crystal Al(lll), Al(llO), and Cu(lll) 



surfaces under surface channeling conditions. The turning point of the trajectory was varied 
by changing the angle of incidence relative to the surface normal. As expected, they find a 
decreasing maximum gy, referred to as "local Fermi momentum", with increasing distance 
from the surface. Their data can be directly compared to our results in sections III Bl and 
111 CI We have displayed the experimental results in the top panels of Fig. |21 for both Al 
surfaces. The contour levels were chosen as to approximately match the experimental data 
point for closest approach to the surface. In the case of Al(lll) the agreement between the- 
ory and experiment is almost perfect whereas results for the Al(llO) surface agree only on 
a qualitative level. This may have two reasons: on the one hand, trajectories along Al(lll) 
surfaces are better defined than along Al(llO). The Al(lll) surface features the "smoother" 
planar surface potential due to its closest packed structure. Additionally, electron spill-out 
and therefore the experimental signal is larger in front of Al(lll) surfaces than in front 
of Al(llO) which may lead to a larger error in the experimentally determined local Fermi 
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FIG. 5: High-momentum tail of p{q\\., z) for z = 2 a.u. The focal Fermi momentum is indicated 
by the dashed line. Contributions above p^p'^ (shaded area) lead to sub-threshold kinetic electron 
emission. 

momentum. 

As the determination of the local Fermi momentum is closely related to measurements 
of "sub-threshold" KE z-dependent momentum distributions are also used for its 

interpretation. In this case, momenta above the local Fermi momentum are related to the 
KE yield (shaded area in Fig. Ej). KE for projectile velocities below (see Sec. |T)) cannot 
be explained within the non-interacting free electron model for metallic conduction bands. 



IV. CONCLUSION 



We have proposed a method to extract information on the surface momentum distribu- 
tions and Compton profiles from ab-initio DFT calculations. Pseudo-wavefunctions from 
such calculations are (partially) Fourier transformed in order to derive the momentum dis- 
tributions p{q) and p{q\\,z) and the quantum mechanical phase-space distributions p^{q,r) 
(Wigner distribution) and p^{q,f) (Husimi distribution). Far from resembling a simple 
Fermi distribution, we find momentum components well above qp originating from correla- 
tion effects and the periodic crystal potential of the solid. Two examples for applications of 



the calculated distributions are presented and compared to experimental results 



2l|. 
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